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ABSTRACT 

The crucial observation on the occurrence of subpulses (overtones) in the Power Spectral Den- 
sity of the August 27 (1998) event from SGR1900+14, as discovered by BeppoSAX (Feroci 
et al. 1999), has received no consistent explanation in the context of the competing theories 
to explain the SGRs phenomenology: the magnetar and accretion-driven models. Based on 
the ultra-relativistic, ultracompact X-ray binary model introduced in the accompanying paper 
(Mosquera Cuesta 2004a), I present here a self-consistent explanation for such an striking fea- 
ture. I suggest that both the fundamental mode and the overtones observed in SGR1900+14 
stem from pulsations of a massive white dwarf (WD). The fundamental mode (and likely 
some of its harmonics) is excited because of the mutual gravitational interaction with its or- 
bital companion (a NS, envisioned here as point mass object) whenever the binary Keplerian 
orbital frequency is a multiple integer number (m) of that mode frequency (Pons et al. 2002). 
Besides, a large part of the powerful irradiation from the fireball-like explosion occurring 
on the NS (after partial accretion of disk material) is absorbed in different regions of the 
star driving the excitation of other multipoles (Podsiadlowski 1991,1995), i.e., the overtones 
(fluid modes of higher frequency). Part of this energy is then reemitted into space from the 
WD surface or atmosphere. This way, the WD lightcurve carries with it the signature of these 
pulsations inasmuch the way as it happens with the Sun pulsations in Helioseismology. It is 
shown that our theoretical prediction on the pulsation spectrum agrees quite well with the 
one found by BeppoSAX (Feroci et al. 1999). A feature confirmed by numerical simulations 
(Montgomery & Winget 2000). 
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1 THE AUGUST 27 (1998) EVENT AND MAGNETAR 
MODEL: CONCORDANCE OR CRISIS 

The lightcurve of the spectacular superoutburst from SGR 1900+14 
in August 27, 1998 exhibited a stable pulsation with period 5.16s 
(Hurley 1999a,b,c; Murakami et al. 1999; Mazets et al. 1999; 
Feroci et al. 1999). Since the modulation frequency is in the 
range of the other three SGRs studied before, Kouveliotou et al. 
(1999) and Hurley et al. (1999a,b,c) concluded that the observa- 
tions provide strong support to the Duncan & Thompson (1992); 
Thompson & Duncan (1995,1996) magnetar model for SGRs. They 
claimed that the observed spindown rate of the pulse period, P — 
1.1 x 10 _10 ss _1 , may be explained by emission of dipolar ra- 
diation from an NS endowed with a very strong magnetic field 
B ~ (2 — 8)10 14 G, a characteristic magnetic field strength in- 
ferred also from the spin down of the pulse period P — 7.47s of 
SGR 1806-20 (Kouveliotou et al. 1998). 

Despite the rough agreement between the SGR1900+14 on 



August 27, 1998 observations and the theoretical prediction of 
the magnetar picture, problems for this model came together with 
that apparent success. On the one hand, it is clear that there is 
an overall consistency between the observations and the magne- 
tar model. However, it is also clear that the simple "giant dipole" 
dynamics is not unique in explaining the spindown history of 
the objects, and in fact several problems with that picture are al- 
ready present when dealing with ordinary, low-field pulsar cousins. 
Among them we can quote the values of the few measured brak- 
ing indexes (which do require modifications from the canonical 
dipole spindown model) and the mismatch between characteristic 
and dynamical-historical ages in the case of some PSR-SN associ- 
ations. Actually, the presence of ordinary pulsars in the high-P, 
high-P "magnetar" region remains puzzling (Manchester 2000), 
and an interpretation of the SGR-AXPs in terms of high (but sub- 
Schwinger) fields is in principle possible (Allen & Horvath 2000). 
On the other hand, it seems that a single population can not account 
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for both ordinary and magnetar objects (Regimbau & de Freitas 
Pacheco 2001). 1 

The event GRB980827 was also detected by BeppoSAX (Fe- 
roci et al. 1999). These observations led to the discovery of an 
extremely regular interpulse set in the X-ray data of GRB980827 
event from SGR 1900+14 (see Fig.l). The interpulses appear sep- 
arated in time ~ 1.1s in between, with no lag. This behavior, Fe- 
roci et al. (1999) advanced, is unexpected and quite difficult to ex- 
plain in the magnetar framework. According to the magnetar model 
for SGRs, global seismic oscillations (Duncan 1998), pure shear 
deformation-induced toroidal modes (standardly labelled as ; T n ) 
with no radial components, i.e., n = overtones, are expected 
to be produced in association with the onset of a new recurrence 
of a "soft" 7-ray repeater. According to Duncan (1998): "... these 
toroidal modes are easy to excite via starquakes because the restor- 
ing force is determined uniquely by the weak Coulomb forces of 
the crustal ions. However, overtones are not allowed because that 
would require far too much energy so as to allow for the extremely 
short period (^ 1 ms, or lower) NS oscillations to be excited". An 
energy that the crust cracking mechanism cannot provide (de Fre- 
itas Pacheco 1998). 

In overall, the simplest and neatest interpretation of that ob- 
servations is that the subpulses in the lightcurve (Power Spectral 
Density) of the burst from GRB980827 are overtones of the funda- 
mental frequency fa = 0.194 Hz (Feroci et al. 1999). As shown 
below, we may be actually detecting the whole WD pulsation spec- 
trum, up to 19 harmonics (Feroci et al. 1999). The importance of 
having a large part of the pulsational spectrum of the object can 
not be overstated. The main purpose of this Letter is to address this 
issue. We show that a self-consistent interpretation for the funda- 
mental mode of pulsation and the harmonics discovered by Feroci 
et al. (1999) could be constructed invoking the excitation of WD 
pulsational p — modes driven by tidal-heating (Pons et al. 2002) 
plus the irradiation (Podsiadlowski 1991) triggered by supercrit- 
ical accretion onto its orbital NS partner in an X-rays ultracom- 
pact ultra-relativistic binary (see a more detailed discussion on this 
model in the accompanying paper by Mosquera Cuesta 2004a). 



2 THE WD EXCITATION ENERGY SOURCE: TIDAL 
HEATING PLUS 7-IRRADIATION 

The basis of the ultra-relativistic compact binary model for SGRs 
(Mosquera Cuesta 2004a), is that during rather sparse catastrophic 
epochs (ATsgrs < lOyr) the WD starts to transfer mass onto a 
low-magnetized (B ~ 10 10 G) rapidly rotating millisecond (ms) 
massive NS (~ 2M0). This process develops via the formation 
of a thick dense massive accretion disk (TDD) very close to the 
innermost stable circular orbit around the NS. The disk becomes 
unstable due to gravitational runaway or Jeans instability, partially 
slumps and inspirals onto the NS. The abrupt supercritical mass ac- 
cretion onto the NS releases a quasi-thermal powerful 7-ray burst 
(GRB), a fireball to say. A parcel of the accretion energy illumi- 
nates with hard radiation (7 — rays) the WD, additionally perturb- 
ing its hydrostatic equilibrium. The WD absorbs this huge energy at 
its interior and atmosphere. 2 This irradiation excites other p-modes 

1 The author truly thanks Prof. J. Horvath (IAG-USP/Brazil) for enlight- 
ning discussions on these issues. 

2 Podsiadlowski (1991; and references therein) has discussed this process 
for irradiated main sequence and evolved stars, where the companion ex- 
pands to a new state of thermal equilibirum, which provides a new mech- 



Table 1. First few radial (1=0) and nonradial (1==0) modes for a 1.05 
M WD model (Temperature =12000 K). Columns 1, 2, 4 and 5, repre- 
sent the 1-value, radial overtone value n, period (s) and frequency (Hz), 
as computed by Montgomery & Winget (1999). 



In P [s] Freq. [Hz] 







1 
1 
1 
1 

2 
2 
2 
2 
2 
3 
3 
3 
3 
3 



of the WD oscillation spectrum, the overtones, since it already pul- 
sates at its fundamental mode because of the tidal interaction of the 
binary (see Table I). 

The abrupt supercritical accretion also perturbs the NS hy- 
drostatic equilibrium which drives it into non-radial nonaxisym- 
metric oscillations that produces GWs due to excitation of the 
NS fluid modes (see Mosquera Cuesta et al. 1998). Further, the 
remaining part of the TDD might be the neighbour environment 
(~ 100km from the NS) where matter carried by the fireball can 
nucleosinthesize to produce the noticeable iron Fe 56 line discov- 
ered by Strohmaier & Ibrahim (2000) during the GRB980827 giant 
outburst. This possibility is explored in a forthcoming communi- 
cation (Mosquera Cuesta, Duarte & de Freitas Pacheco, in prepa- 
ration). (We address the reader to the related paper by Mosquera 
Cuesta 2004a, for a full description of the interacting relativistic 
compact binary here pictured). 

2.1 A possible origin for the fundamental mode frequency 
and overtones 

In this section we argue that the theoretical modeling of WD pul- 
sations as performed by Montgomery & Winget (1999) can accu- 
rately account for the spectrum of frequencies discovered by Bep- 
poSAX during the event GRB980827 from SGR1900+14 (Feroci 
et al. 1999). We note in passing that although the Montgomery 
& Winget (1999) numerical simulations focussed on the study of 
long period gravity-induced WD oscillations, which are highly rel- 
evant for the potential identification by the Whole Earth Telescope 

anism to drive mass transfer onto the NS. This alters the binary evolution, 
and may bring a new evolutionary stage during which the orbital period in- 
creases, leading to larger orbital period during and at the end of the mass 
transfer. Although the WD in this relativistic binary is a degenerate com- 
pact star with an atmosphere, we believe a similar behavior should also take 
place in it, since the new thermal timescale, <x Ai? WD , is comparable to 
the heat (7-rays) difussion time in the outer layers. 
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Table 2. A subset of the WD frequencies displayed in Table I, versus 
the modulation spectrum from SGR 1900+14 on August 27 (1998) as 
discovered by Feroci et al. (1999). 



Overtones. 


Num. Model [Hz] 


BeppoSAX [Hz] 


Mismatch (%) 


1= 0, n = 1 


0.1846 


0.194 


5.0 


1= l,n= 1 


0.3898 


0.389 


0.2 


1= 2,0, n = 2,3 


0.7644, 0.7845 


0.775 


1.0, 1.23 


1= 1,2, n = 3,3 


0.9025, 1.0068 


0.969 


4.54, 3.78 


1=1, n = 4 


1.1430 


1.161 


1.55 



of pulsating ZZ-Ceti WDs in surveys, those models also provide the 
spectrum of pulsations of the short period pressure-driven (p-mode) 
pulsations. 

Any perturbation of the hydrostatic equilibrium of a canonical 
WD will grow on its dynamical timescale 
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The WD normal mode spectrum (p — modes) is obtained 
from the radial wavenumber k r defined by (Montgomery & Winget 
1999) 



kr = (l/o" 2 Cs)(o- 2 - L 2 )(a 2 



N 2 ) 



(2) 



with a the mode angular frequency and c s the sound speed 
in the star material. Here the squared Lamb/acoustic frequency is 
defined 



1(1 + 1) 



(3) 



where r is the radial variable, and A^ 2 the Brunt-Vaisala fre- 
quency. For p — modes: 

2 



2 . r 2 



N 



and 

k r ir 



r 2 dr ■ 



(4) 



(5) 



where ri, n are the inner and outer turning points, respec- 
tively, at which k r = for a given a. 

Montgomery & Winget (1999) studied the pulsational modes 
of a massive WD (Mwd ~ I.IMq) which possess a hydrogen 
atmosphere of about M^rf 4 ~ 1O" 6 M , with and without core- 
crystallization. The full spectrum of one of their numerical models 
is displayed in Table I, and a selected subset is displayed in Table 
II 3 to compare with BeppoSAX observations (Feroci et al. 1999). 



3 DISCUSSIONS 

3.1 The WD required mass 

If we use the picture being introduced here to explain the pulsa- 
tion discovered in GRB980827 (Hurley 1999a; Kouveliotou 1999a; 
Murakami et al. 1999; Feroci et al. 1999), it is easy to see that the 
WD mass (see Table 3) needed to produce pulsations with timescale 

3 A curious note on this result is that it exhibits the same number of over- 
tones as the one Feroci et al. (1999) found in the GRB980827 data from 
SGR1900+14. 
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Figure 1. Graphic comparison between the observed modulation (main six 
overtones) in SGR1900+14 as detected by BeppoSAX (Feroci et al. 1999) 
and the theoretical pulsation spectrum (vertical dashed lines) computed by 
Montgomery & Winget (1999). We claim that the overtone of frequency 
near 0.96 Hz might have been enhanced in power via irradiation during the 
7-rays superoutburst. 



of 5.16s is ~ 1.1 M (Montgomery & Winget 1999; Shapiro & 
Teukolsky 1983). The interpulses discovered by the Italian team 
of the BeppoSAX collaboration can also be explained in a simple 
manner (except that there are instrumental errors): they are a kind of 
WD ringing overtones. The full attenuation of the oscillations will 
be achieved on a much longer timescale via the wave-leakage of 
radial modes as discussed by Hansen, Winget & Kawaler (1985), 
see below. So the observed pulsations may stem from a WD, not 
from a slowly rotating hypermagnetised NS. 



3.2 SGRs spindown rate 

Moreover, the SGR 1900+14 spindown determined by Kouve- 
lioutou et al. (1999) and the pulse period increase found by Mu- 
rakami et al. (1999); Marsden, Rothschild & Lingenfelter (1999); 
Harding, Contopoulos & Kazanas (1999) may be explained as fol- 
lows: because the WD is a gravothermal system, as soon as it loses 
mass (when overflowing its Roche lobe) its negative specific heat 
forces it to expand until a new dynamical equilibrium radius is 
found. Then, the next stage of pulsation should occur with a slightly 
longer period compared to the previous one (secular increase). As 
a result, the pulse period increases. However, the interplay between 
the WD thermal cooling during the post-outburst reemission (which 
implies post-burst shrinking) and tidal heating via gravitational in- 
teraction with the NS (which implies pre-burst expansion) may lead 
to a slight up-down-up change in the fundamental mode of pul- 
sation. This processes may help to understand the observed tiny 
changes in the SGR 1900+14 period: near superoutburst it is a bit 
longer, while in quiescence it shortens. 

Thus, the rate of variation of the viscously attenuated pulsa- 
tion time, T pu i se — 5.16 s, divided by the timescale t w -i ~ 10 3 yr 
needed to dissipate via wave-leakage for radial modes 4 the ab- 



4 Physical support for the figures that have been used here comes from 
the detailed numerical calculations of pulsation periods for DA and DB 
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Table 3. SGRs observed pulsation periods and in- 
ferred masses for the WD in the binary model. 



SGR 


P[s] 


Mass [Mq] 


1900+14 


5.16 


1.1 


0526-66 


8.1 


0.70 


1806-20 


7.47 


0.80 


1627-41 


6.7 


0.95 



sorbed energy ~ 10 erg, gives us the observed "mean" spindown 
rate in SGR 1900+14 (here we assume this the dominant source 
available to the WD) 

/ dr^\ 5.16 s = 10 ! 

V dr w -i ) 3.17 x 10 10 s 

In our view this is the origin of those important time varia- 
tions. Discrepancies might be caused by the differences in the WD 
model used (see Table 1). We stress that in a realistic star the over- 
tone frequencies depend on the WD matter EOS and the degree of 
core crystallization. Changes in pulsation (fundamental) periods of 
~ a few percent are expected to occur once lattice crystallization 
onsets while the harmonics (overtones) remain almost unchanged 
(Montgomery & Winget 1999). 



3.3 X-ray emission during quiescent, pre- and post-burst 

What about X-ray emission and pulsation in long-term quiescence, 
pre- and post-burst? A number of viable processes may explain 
why SGRs glow in hard X-rays over some months before and af- 
ter undergoing dramatic transients such as GRB980827 in SGR 
1900+14, as well as the long-term quiescent emission: 1) changes 
in the gravitational potential at each orbital revolution (Podsiad- 
lowski 1991,1995), 2) crust cracking driven by the WD strong 
magnetic field (analogously to the NS case, see de Freitas Pacheco 
1998), and 3) plate tectonics (Rothschild, Marsden & Lingenfelter 
2001) induced by any or both the mechanisms just listed. 

In this paper, as a first approach, we shall focus our discus- 
sion on this issue in the context of the scenario number 2) above. 
The remaining possibilities will be addressed elsewhere. This key 
property can be explained if we suppose that the quiescent soft X- 
rays luminosity Lx GR " ~ 10 35 erg s" 1 (Kouveliotou et al. 1999) 
is powered by the release of WD crustal elastic energy , inasmuch 
as in the current picture of magnetars (de Freitas Pacheco 1998). In 
this case, we get 

ga ; r/J?c ° re x 4irR 2 WD x AR WD > L s x GRs x r (7) 

where r ~ 10 4 yr is the system lifetime since its formation, 
B sur f ~ 10 9 G WD crustal magnetic field, as in the super magne- 
tised WD: RE J0317-853, Mq = 1.35, P = 725 s (Wickramasinghe 
& Ferrario 2000), and B core ~ 10 13 G is the WD core magnetic 



WDs given by Hansen, Winget & Kawaler (1985), where timescales for 
dissipation of energies such as the one absorbed by the WD in this scenario 
are estimated for wave-leakage via unstable radial modes with e-folding 
timescale m ~ 10 3 yr, for nonadiabatic driving of the mode with energy 
10 44 erg. 

5 An original suggestion by Prof. Posiadlowski (2000), Dept. of Physics, 
Oxford University, England, UK. 



field. This yields a WD crust thickness: ARwd ~ 70 km 6 . The 
crust-cracking induces shear stresses which dissipate energy caus- 
ing excitation of WD oscillation modes quite similiar to the NS 
case (de Freitas Pacheco 1998, and references therein) 7 . We thus 
emphasize that the conclusion by Mazets et al. (1999): "... the pro- 
cesses accounting for emission of the narrow initial pulse and the 
long pulsating tail in both SGR 0526-66 and 1900+14 are sepa- 
rated in the source not only on time but in space..." , is realized in 
the picture introduced here, i. e., the NS releases the superoutburst 
while the WD the subsequent tail of pulsations via the mechanism 
invoked above. 



4 CONCLUSIONS 

As a summary, the NS low magnetization in this model leads 
SGR1900+14 (and perhaps all SGRs) under the pulsar death line 
making it undetectable as a binary radio pulsar, a viewpoint con- 
firmed by Xilouris et al. (1998). Overall, since several SGRs are 
enshrouded by intervening galactic dust and gaseous nebulae (SNe 
remnants: Gaensler et al. 2001), optical observations of the sug- 
gested hot WD (~ 12000 K) may render a breakthrough. This 
fact makes it a systematic high resolution search for such opti- 
cal (or infra-red) counterparts of SGRs a timely and promissing 
task. In this line, we remark that Ackerlof et al. (2000) have pur- 
sued optical follow-ups of SGRs over 10 outbursts: 8 events from 
SGR1900+14 and 2 events from SGR 1806-20. Although a care- 
ful search for new or variable sources in the view field of those 
SGRs was performed, no optical counterparts were seen down to 
mv ~ 15. Thus the search for such SGRs counterparts remains, 
and perhaps VLA, GEMINI or KECK observations in the K-band 
could unravel them. On the other hand, it will be extremely rele- 
vant if the propeller model of Marsden, Rothschild & Lingenfelter 
(2001) and Alpar (1999, 2000) could provide a clean explanation 
for the observed spectrum of pulsational modes. This may help to 
discriminate among the SGRs accretion models through future ob- 
servations. 

We deeply thank Mike Montgomery and Don Winget for 
kindly running their numerical code to simulate WD non-radial 
p-mode pulsations on our request, and also for the many fruitful 
discussion on this subject. 
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